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Abstract. We discuss the possibility of simultaneous and sequential synchronisation in vertical and hor- 
izontal arrays of unidirectionally coupled discrete systems. This is realized for the specific case of two 
dimensional Gumowski-Mira maps. The synchronised state can be periodic, thereby bringing in control 
of chaos, or chaotic for carefully chosen parameters of the participating units. The synchronised chaotic 
state is further characterised using variation of the time of synchronisation with coupling coefficient, size 
of the array etc. In the case of the horizontal array, the total time of synchronisation can be controlled by 
increasing the coupling coefficient step wise in small bunch of units. 

PACS. 05.45.Xt Synchronization; coupled oscillators - 05. 45. -a Nonlinear dynamics and nonlinear dy- 
namical systems 



■ 1 Introduction nised behaviour of neuronal networks [SJE]. Moreover pat- 

~,] . terns of synchrony and phase shifts in animal locomotion 

a- 

Synchronisation of the dynamical variables of coupled sys- } s gaining importance as a field of active study 



m a mi mi 



terns is an important nonlinear phenomenon where intense |TT1[T2] . In general, the synchronised networks for analysing 

research is being concentrated recently 0. This is prob- or modelling all these physical or biological situations arc 

ably because of its engineering applications like spread constructed by coupling basic dynamic units with a well 

spectrum and secure data transmission using chaotic sig- defined connection topology that can be nearest neigh- 

nals HEE], control of microwave electronic devices 4 , bour, small world, random or hierarchical architectures, 

graph colouring etc. Also communication between differ- T n addition, in specific applications like communication or 

ent regions of the brain depends heavily on the synchro- neural networks, a realistic modelling may require the in- 
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transmission or processing speeds among the units ,13- In etc. These two arrays can be further worked together to 

any case, it is found that the collective dynamics depends produce square lattice networks with desirable or useful 

crucially on the connection topology |14| . inter connections. 

The simplest yet the most widely used topology in this The array is realised here with a two dimensional dis- 

context is the linear array and its combinations. The study Crete system or map as the local unit and a connection 

of synchronisation in arrays of systems was first applied to that involves a non linear function forming part of the map 

laser systems |15U16 | which has relevance in optical com- function. The stability of the simultaneously synchronised 

munication systems. Since then the occurrence of synchro- state for the vertical array is studied by computing the 

nisation in coupled map lattices has been extensively stud- Maximum Conditional Lyapunov Exponent (MCLE) |19| . 

ied with many consequent applications |17) . Such systems, so that the minimum coupling coefficient required for on- 

with synchronisation in temporally chaotic but spatially set of synchronisation can be deduced. The dependence of 

ordered units forming an array, is applied in many situ- the characteristic response time r s on the coupling coef- 

ations like data driven parallel computing |18| . However ficient e is analysed numerically. A horizontal array with 

most of these cases studied so far involve continuous sys- the same dynamics is constructed with each unit driven 

terns of chaotic oscillators. by the previous one, modelling an open flow system and 

T ,, . . , , , leading to sequential synchronisation. In this case the time 

In this paper, we consider two such regular arrays, one 01 

, . , 1 ,, ,•, , . ,1,1, 1 1 , 1 taken for the last unit to synchronise is taken as the total 
vertical and the other horizontal, that works under the 

, . , , ,, .,. response time r,. The behaviour of its average for differ- 

drive response mechanism, where the connection is umdi- r 

1 m c j ii i ii c i j i u ent initial conditions and size N of the system are studied, 

rectional. We find that the former setup leads to simulta- J 

i ... , ., ,, , , , ,. . The additional time or delay time 77 required for the last 

neous synchronisation while the latter results 111 scquen- J 

,. , , ... TT iii-i ,, unit to synchronise after its previous one has synchronised 

tial synchronisation. Here we would like to comment that J 

n , i , j . ! , , , is found to saturate with system size. Moreover we note 

m most ot the connected networks, the synchronisation J 

. r , , i tt ^1 1 an interesting bunching effect where the total t, can be 

is found to occur simultaneously. However the topology 00 

, 1 i • , . , , . , , i , , , controlled by varying the value of e in bunch of m units, 

m the linear horizontal array introduced here develops 

synchronisation sequentially and the delay time from one In Section 2, we introduce the basic unit which serves 

unit to the next can be adjusted by external control. This as the driving as well as the driven systems with identical 

mechanism therefore would be useful for many technolog- individual dynamics. The concept of generalised synchro- 

ical applications. These two types of synchronisations are nisation and its stability in the context of unidirection- 

characterised using response time (which is the time for ally coupled systems is also discussed. The construction 

synchronisations to stabilise), size effect, bunching effect and the collective dynamics of the vertical array and the 
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characterisation of simultaneous synchronisation is given or horizontal array and driven unidirectionally, they are 
in section 3. In section 4, we introduce sequential syn- capable of synchronising to the same chaotic state, 
chronisation and its control due to the bunching effect In the context of unidirectionally coupled systems, the 
of the unidirectionally coupled units. Our concluding re- type of synchronised behaviour called generalised synchro- 
marks are given in section 5. nisation has been attracting much attention recently [231 

124) . Here the states of the driving system Xd and the 



2 Basic Dynamical unit and Generalised 
Synchronisation 

The basic unit used here for the present analysis of syn- 
chronisation in arrays is a two dimensional discrete sys- 
tems, which serves both as the driving and driven systems 
defined in the phase space X(n) — (X(n), Y(n)). The spe- 
cific system chosen for this work is the Gumowski-Mira 
recurrence relation |2U| given as 

X(n + 1) = Y(n) + a (l - 6Y(n) 2 ) Y(n) + f{X(n)). 



driven system Xdr are dynamically related by a function 
F such that the relation Xdr(t) = F(Xa(t)) is true once 
the transients are over. The form of F can be smooth or 
fractal and in either case, the procedure for finding the 
same can be complicated. Hence often an auxiliary sys- 
tem identical to the driven system is introduced as X a (t). 
The initial conditions of Xd r and X a are taken different 
(both being individually chaotic in dynamics) but lying in 
the basin of the same attractor. Once the transients have 
settled, the dynamical equivalence of Xd r (t) and X a (t) is 
taken as an indication of generalised synchronisation be- 



Y(n + 1) = -X(n) + f(X(n + 1)). (1) tween X d {t) and X dr (t). 



2{1 - ii)X 2 (n) 
1 + X 2 (n) 
n refers to the discrete time index 



where f(X(n)) = fiX(n) H ■ v2 , x — and 



3 Simultaneous Synchronisation in a Vertical 

Our earlier investigations in this system reveal that Array 



<|TTp is capable of giving rise to many interesting two di- 
mensional patterns in (X, Y) plane that depend very sen- 
sitively on the control parameter \i |21|. This can be ex- 
ploited in decision making algorithms and control tech- 
niques for computing and communications. We have tried 
three different coupling schemes in two such systems |22| 
and found that they are capable of total or lag synchroni- 
sation in periodic, quasi periodic or chaotic states, when 
N such systems are geometrically set to form a vertical 



We extend the above concept to construct a vertical array 
of N identical systems, [X^ r (n), X^ r (n) ■ ■ ■ X d w r (n)] ; each 
driven independently by Xd(n). All the systems are iden- 
tical and individually evolve according to Fig. show 
the above scheme of construction of vertical arrays. 
Here the driving system follows the dynamics 

X d (n + 1) = Y d (n) + a (l - bY^n)) Y d (n) + f(X d (n)). 
Yd{n + l) = -X d (n) + f(Xd(n + l)). (2) 



4 
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dr 



Fig. 1. Schematic view of the construction of a vertical array 



with f{X d (n)) = fi d X d {n) + 



2(1 - fi)Xj(n) 
1 + Xj(n) ■ 



The i driven unit in the vertical array has the dynamics 



X dr (n + 1) = Yl(n) +a(l- bYJ r (»)J Y^ r {n) 

+ f(X dr (n)) + e (f(X d (n)) - f(X dr (n))) 
Yl{n + 1) = -X dr (n) + f(X dr (n + 1)). (3) 

with f(X dr (n)) = MrX dr {n) + 2{ \^ d ;f.\ {n) where 
e is the coupling coefficient of the unidirectional coupling 
applied to the X variable through the function f(x). The 
parameters a and b are set as a = 0.008 and b — 0.05. The 
total number of units considered is N = 51. The value of 
fi dr is chosen to be the same for all the 50 driven units. 
We can realise synchronisation for different combinations 
of values of fi dr and fi d with fi d , in general different from 
fi dr . For the special case of fi d = fi dr all the 51 units 
synchronise including the driving system, when started 
with different initial conditions. 

Fixing the value of coupling coefficient e = 0.9, the 
values of ^i dl [i dr for which synchronisation is feasible in 
the 50 driven systems is plotted in fig. [5] In the param- 
eter plane considered here in the range —0.2 < \L d < 
—0.5, —0.2 < /idr < —0.5, the points marked * indicates 
(l^d, A*dr) values leading to synchronised periodic state with 



H □ 1 □ 

□ □ 



-0.5 
-0.5 



Fig. 2. Points in the fid — [idr plane for which synchronisation 
of the driven systems are possible with e = 0.9. Points marked 
* indicates fid, fidr values leading to synchronised states with 
periodicity less than 15. Points marked □ correspond to syn- 
chronisation in higher periodic states or mostly chaotic states. 

periodicity less than 15. Points marked □ indicates syn- 
chronisation in higher periodic state or mostly chaotic 
states. 

For specific cases like fi d = fi dr = —0.39 both the 
driving system and driven systems are in chaotic state in- 
dividually. With e — 1.56 all the 50 driven systems are 
synchronised in the chaotic state while the driving system 
is asynchronous with them. However when e is slightly in- 
creased to 1.6 all the 51 units are found to synchronise in 
the chaotic state. Fig. gives this chaotic synchronisa- 
tion between two participating driven systems for e = 1.6, 
where the iterates of the X variable of the 6 th and 49 th 
units are plotted, after the transients have died out. 

For fi d = fi dr = —0.23, individually the systems are 
chaotic. For e = 0.9, all the N driven systems are syn- 
chronised to the same periodic state of periodicity 15. But 
the driving system is also synchronised only when e is in- 
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For Hd = —0.2, the driving system is in periodic 8 cy- 
cle. For fidr = —0.39 the driven systems are individually 
chaotic. For e = 0.9 all the N driven systems are synchro- 
nised in the chaotic state. Fig. 0] shows the synchronised 
chaotic state for the above case. Here the iterates of the 
X variable of the 48 th unit and 10 th unit are plotted after 
the transients have died out. 




(a) 



49 
X i 



-4 -3 -2 -1 1 2 3 4 5 6 



(b) 

Fig. 3. (a) Chaotic synchronisation between two participating 
driven systems with e — 1.6, fid = fidr = —0.39. Here the 
iterates of the X variable of the 6 th and 49 th units are plotted, 
(b) Synchronised periodic 15 cycle for fid = fidr = —0.23 with 
e = 1.6. The iterates of the X variable of the driving system 
and 49 th driven unit is plotted. 

creased to 1.6. Fig. [^Jd gives this synchronised periodic 15 
cycle for e — 1.6, where the iterates of the driving system 
and the 49 th unit are plotted. Thus the driving system 
and the driven system can be simultaneously synchronised 
only when fid — fidr and when e is very large ie e ~ 1.6. 
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Fig. 4. Synchronised chaotic states between two driven sys- 
tems for fid = —0.2, fidr = —0.39 with e = 0.9. Here the 
iterates of the X variable of the 48 th unit and the 10 th units 
are plotted. 



The condition for the stability of generalised synchro- 
nisation is discussed using the Maximal Conditional Lya- 
punov Exponent Xmcle |25U26| . Here the Lyapunov Ex- 
ponent of the driven system is calculated and it is differ- 
ent from the uncoupled system, since it depends on the 
dynamics of the driving system also. The condition for 
the stability of generalised synchronisation is that Xmcle 
should be negative |27|. From equations J2Jl and @ the 
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Jacobian matrix for the i th unit can be written as 



M 



V 1 L )dX' 8Y* 



f)G' 
8X ' 



dG' 



(4) 



For the above case the coupling coefficient is varied 
in steps from 0.84 to 0.97 and the time taken for reaching 
synchronisation in the driven systems is noted. Fig. Ogives 



where X l dr (n + 1) = F^X, Y) and Yj r (n+1) = G^X, Y). the variation of thus average response time t s (averaged 
If a 1 and a 1 are the eigen values of the product of the 



Jacobian matrices at every iteration such that a 1 > a 2 , 
then [2HI 

Amclb = lim — ln\a l \ (5) 

^mcle can be calculated numerically for different e values 
using JSJl. 

We consider the case Hd — —0.2 and /idr — —0.39 
which we have discussed above and calculate \mcle for 
different e values. Calculations are done for 10000 iterates 
after leaving initial 70000 iterates as transients. In fig. [5] 
the values of \mcle for different e values are plotted. It 
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Fig. 5. The variation of Maximum Conditional Lyapunov Ex- 
ponent (Amclb) of the driven system with the coupling coeffi- 
cient e. Here fj-d = —0.2, Hdr = —0.39. The minimum coupling 
coefficient for synchronisation e m i n = 0.829. 



0.82 0.84 0.86 0.88 0.9 0.92 0.94 0.96 0.98 

£ 

Fig. 6. The variation of the response time (r s ) which is the 
total time for synchronisation with coupling coefficient e for 
the 50 th unit. (r s ) is almost constant for values of e > 0.87. 

over 10 different initial conditions) with e for the 50 th unit. 
It is interesting to note that (r s ) is almost constant for 
values of e > 0.87. In this case since the synchronisation 
is simultaneous and coupling is unidirectional and similar, 
the average (t s ) is independent of the size of the array N. 

4 Sequential synchronisation in a Horizontal 
Array 



In this section a horizontal array of N identical systems 
with open ends, where each unit is driven by the previous 
one is introduced. The coupling is through the nonlinear 
is found that, Xmcle crosses zero at e = 0.829, which is function f(X,Y) as in the previous case. Fig. [7|gives the 



the minimum value of e |19j viz e m j n such that for e > e n 
the synchronised state is stable. 



schematic view of unidirectional coupling in a flow which 
consists of N units. 
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SI 



S2 



• S3 



Fig. 7. Schematic view of unidirectional coupling in a flow 
system. 

The i th unit in the horizontal array follows the dynam- 
ics 

X\n + 1) = Y l (n) +a(l- bY 1 * (n)) Y^n) 

+ f(X i (n)) + e{f(X i - 1 (n))-f(X i (n))) 
Y\n + l) = -X l {n) + /(A> + 1)). (6) 

with f(XHn)) = ftXUn) + 2(1 ~ M) f ' S n) ■ 

The control parameter zx in the same for all the units 
such that the units are chaotic individually. This set up 
is found to give rise to sequential synchronisation in the 
array. In our calculations we consider an array of 51 units. 
This can be extended to any number of units N. 

For /x = —0.23 where the individual systems are chaotic, 
and coupling coefficient e = 1.9, we find that synchroni- 
sation sets in sequentially with the 2 nd synchronising af- 
ter the first, the third after the second and so on. The 
time taken by the last unit to synchronise is taken as (r s ) 
which is the average total response time for the whole ar- 
ray. Fig. [5] shows this synchronised chaotic state after the 
last unit has synchronised. The (r s ) is found to vary with 
coupling coefficient e as shown in Fig. [5| It is found that 
(t s ) has a minimum value for a particular e which in this 
case is e = 2. 

The delay time 77 ie, the additional taken for the iV th 
unit to synchronise after its previous one has synchronised 
is defined as 77 = rf — rf -1 . This t; is found to saturate 




X 1 

Fig. 8. Synchronised chaotic state for [i = —0.23 with e = 1.9. 
The iterates of the X variable of the first unit and the last unit 
(51 st unit) are plotted. 




1.9 1.92 1.94 1.96 1.98 



2.02 2.04 



Fig. 9. The variation of total response time (t s ) of the 51 
unit with coupling coefficient e. (r s ) has a minimum value for 
e = 2. 

with the system size as shown in fig.EH Beyond N = 
35, 77 is almost constant. 

An interesting observation in this horizontal array of 
units is a bunching effect that reflects in the total response 
time (t s ). For this instead of fixing the same value for the 
coupling coefficient e for all the units, we fix its value for a 
particular number of units and increase it in steps for the 
next bunch and so on. Then the total (t s ) is found to be 
smaller compared to the previous case of the same e for 
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Fig. 10. The delay time n which is the additional time taken 
for the N th unit to synchronise after its previous one has syn- 
chronised is found to saturate with system size N. Beyond 
N = 35, Ti is almost a constant. 



Fig. 11. Change in the response time (r s ) of the last unit with 
bunch size m. Here fi = —0.23. Curve a is for e m ax = 1.91, when 
the bunch size m = 7, the response time (r s ) is minimum, b is 
for £ m ax = 2.01, the response time (t s ) is minimum for m = 8. 
all the units. Moreover this time depends on the size of In both cases (t s ) is found to be less than the case when we 
the bunch and is minimum for a certain number of units apply e max to all the units, 
in each bunch. 

We report a few specific cases. With \x = —0.23 the 
value of e is increased in steps of 0.001 for each bunch 
so that e for the last bunch is e max = 2.01 for different 
bunch sizes. In each case the total response time (r s ) is 
found. Fig. 1111 shows how the response time (r s ) changes 
with the variation in the bunch size to, ie., number of 
units in each bunch. The response time (r s ) is minimum 
when the bunch size to = 8. For m = 8, (t s ) — 138612 
iterations, whereas when e = e max = 2.01 for all the units, 
the response time (r s ) = 144404 iterations. 

As a second case for same /i = — 0.23e max is taken as 5 Conclusion 
1.91 and calculations repeated as above. In this case (r s ) 

is found to be minimum and is 152150 iterations when In this work we report how synchronisation in an array 
the size of the bunch is m — 7 as shown in Fig. ^2 If of systems can be made more efficient and flexible to suit 
e = 1.91 for all the units, (t s ) is 161953 iterations. We specific applications. We consider two such arrays, vertical 



to bunching must be reflected in the response time of each 
bunch. So for the minimum case, the response times for 
the last unit of the first bunch (ie., 8 th unit), last unit 
of the second bunch (ie., 16 th unit) etc. are noted with 
bunching. The same quantity with e same for all the units 
ie., without bunching are also noted. (t s ) thus obtained 
are plotted against the respective units in Fig. It is 
found that except for the 8 th unit in the first bunch the 
response time is less in the case of bunching. 



observe that the decrease in (r s ) for the whole array due and horizontal, working under the drive-response mecha- 
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Fig. 12. The response time (r a ) is plotted for different units 
for n = -0.23, £max = 2.01. The dotted line gives the total to Saturate with S y stem size 



to overcome the transients is found to saturate beyond 
a certain value of e. The horizontal array exhibits many 
interesting features useful for technological applications. 
In this case the synchronisation sets in sequentially from 
unit to unit along the array since the coupling is unidirec- 
tional. The total response time for the whole array has a 
minimum as e is varied. The additional time required for 
the last unit to synchronise after the previous one is found 



response time (r 3 ) of 8 th unit, 16 th unit, 24 th unit etc. when 
£max = 2.01 is applied to all the units. Full line gives the total 
response time (r 3 ) of 8 th unit (last unit of 1 st bunch), 16 th unit 
(last unit of 2 nd bunch), 24 th unit (last unit of 3 rd bunch) etc. 
when e is increased in steps for each bunch so that e max = 2.01 
for the last bunch. Here 2.005 < e < 2.01 with step size 0.001. 
Thus bunching can control the total response time (r s ) of the 
array. 



nism with a two-dimensional discrete systems as the unit 
dynamics. 

It is observed that synchronisation sets in all the sys- 
tems simultaneously in the vertical setup. The minimum 
value of the coupling coefficient e required for stability of 
synchronisation is computed numerically from the Maxi- 



We further note that the total response time for the 
whole array can be reduced by introducing bunching with 
step wise increase of e from bunch to bunch. There exists a 
specific bunch size giving minimum time which depends on 
the choice of the parameter and the maximum e given to 
the last bunch. This makes the sequential synchronisation 
flexible and controllable to suit specific applications. 

At present we do not find any specific reasons for the 
above findings to depend on the unit chosen. For differ- 
ent choices of the unit dynamics the behaviour should be 
qualitatively similar. To establish this generality and ap- 
plicability of the present technique, we are trying it out 
for a number of systems. The results will be published 
elsewhere. 



mum Conditional Lyapunov Exponent. The possible choices 
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